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Abstract 

In this letter we introduce a possible measure of the size of CP violation in the Stan- 
dard Model and its extensions, based on quantities invariant under the change of weak 
quark basis. We also introduce a measure of the "average size" of CP violation in a 
model, which can be used to compare the size of CP violation in models involving extra 
sequential or vector-like quarks, or left-right symmetry. 

PACS: ll.30.Er, 12.15.Ff, 12.60.-i, 14.65.-q 

The definition of "Maximal CP Violation" and the related problem of finding an adequate 
measure of the size of CP violation in the Standard Model (SM) has been object of interest 
j|, ^, ||, |], H since experiment || revealed that the phase appearing in the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix had to be large to explain the observed CP violation in the K°- 
K° system. As shown by Wolfenstein Q any definition of maximal CP violation based on 
maximizing the phase appearing in the CKM matrix does not make sense because it depends 
on the parametrization of the CKM matrix itself. A phase of 7r/2 in one parametrization does 
not correspond to a phase of tt/2 in another parametrization. Gronau and Schechter [0] used 
the construction of the CKM matrix as a product of three 2x2 unitary matrices (i. e. the 
Murnaghan construction) to find a certain combination of phases (the "invariant phase" ) that 
remained invariant under most rephasings of the quark fields. However, this invariant phase 
depends on the adoption of the Murnaghan construction as well as on the order in which the 
2x2 matrices are multiplied. The formulation in terms of rephasing invariants of the CKM 
matrix was carried out independently by Jarlskog ||] and Dunietz, Greenberg and Wu H 
leading to two different definitions of measures. In Refs. Q ||] the condition of maximality is 
that some of the quantities acp = 21m (a/3*)/(|a| 2 + |/3| 2 ), with either a = VijVki, (3 = VkjVu 
or a = VijVft, j3 = Vp-jV^ or a = VijVj^, = VuV^ are maximal, Vij being the ij element 
of the CKM matrix. Note that the numerators of all the quantities acp are equal up to an 
overall sign by the unitarity of the 3x3 CKM matrix but the normalization differs. In Ref. || 
CP violation is said to be maximal when the products Im VijV^VkiV^ acquire its maximum 
absolute value. This occurs when there is maximum mixing in the CKM matrix: the modulus 
of all the matrix elements is l/v3 and the phase of V u dV* d V cs V* s is 27r/3. Nevertheless as 
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pointed out by Botella and Chau Q this simple and elegant definition cannot be generalized 
to a higher number of quark generations because unitarity of a N x N CKM matrix does not 
imply that all the quantities Im VijV^VkiV^ are equal in modulus, and one has the case that 
when maximizing two of them a third one becomes zero. In this letter we follow the work of 
Refs. |3], ||, J5| although our definitions will be based on quantities invariant not only under 
rephasings of the CKM matrix but under arbitrary quark basis transformations. 

The formulation of CP violation in terms of quantities invariant under a change of weak 
quark basis is complementary to the usual formulation in terms of CKM phases, and seems 
more adequate to define a measure of CP violation. In Ref. [|]] Jarlskog realized that, if 
properly normalized, the quantity Im det [M u M\, M^mJ], invariant under weak quark basis 
transformations, can be looked as a measure of CP violation in the SM. (M u ^ are the up and 
down quark mass matrices respectively.) All CP violation observables depend in a perhaps 
complicated way on this quantity, which then gives us a measure of CP violation in the SM. 
However, the normalization is not uniquely determined and can lead to different measures. 

To extend this formulation to a more general model, let {M{\ be the set of quark mass 
matrices of the model, together with their transpose, conjugate and hermitian conjugate. 
These matrices are such that under a change of weak quark basis Mj — > uj Mill- with f/j, U- 
unitary matrices. For instance, in the SM it is enough for our purposes to choose {Mi} = 
{M u , Mdi M|, Aft}. We can construct a quantity invariant under a change of weak quark 
basis by taking the trace of a product of n matrices, tr Mi ■ ■ ■ M n , which can be repeated, 
chosen in such a way that Ui+i = U-, Vi (with U n +i = U\). The number n of matrices may 
not be arbitrary in specific models. Then CP conservation (for the detailed proof see Ref. 
|J) implies that the imaginary part of this trace, which we will call an "invariant" from now 
on, is real: 

I = Im tr Mi • • • M n = . (1) 

A set {/} of invariants can be constructed in this way. The vanishing of all these invariants is 
a necessary condition for CP conservation. Furthermore, in specific models sets of invariants 
can be found with the property that their vanishing is also a sufficient condition for CP 



conservation |4|, 1C, 11, 12]. For instance, in the SM the condition Im tr H^H^HuH^ = 
with H x = M x Mlis sufficient for CP conservation [10] (note that in the SM Imdet[H u , H<j\ = 
21m tr H^HdHuHj ). In the SM with an arbitrary number of generations of quarks all the 
invariants are constructed with products of the hermitian matrices H u ^- Any CP violation 
observable must depend, perhaps in a complicated way, on these invariants. Generalizing 
the results in Refs. [||, we argue that these invariants properly normalized ("reduced 
invariants") give measures of CP violation. 

As pointed above, the normalization is not unique. Nevertheless, by imposing some 
properties on the reduced invariants we can determine a plausible form. Let us consider the 



2 



invariants 



h 
h 



Im tr H 2 u H d H u H 2 d , 
lmtiH^H d H u Hj. 



(2) 



In the SM, all the invariants are proportional to I\, which is the lowest order nontrivial one. 
Explicit calculation shows that I2 = (m^ + m 2 + mf) Ii, with rrii the mass of the quark i. If 
we require that the reduced invariants I±, I2 are equal, the relative normalization must be 
( m u + m c + m t ) = tr H u . The same can of course be done with 73 = Im tr H 2 H d H u H^. 
This suggests that the adequate normalization for an invariant / with n u powers of H u and 
n d powers of H d should be I = //(tr H u ) nu (tr H d ) nd . This generalizes to other models as 
follows: for an invariant I = Im tr M\ ■ ■ ■ M n the reduced invariant / is 



This normalization has the remarkable property that — 1 < I < 1 for any reduced invariant I 
and arbitrary matrices M\, . . . ,M n , as holds for the CP asymmetries acp discussed in Refs. 
[Q, ||] . This property can be easily shown using the Schwarz inequality for the scalar product 
of matrices (A, B) = tr AB^ which in this case reads (tr AB 1 ") 2 < tr AA^ tr BBK 

A deeper reason to choose this normalization is that for a measure of CP violation we 
want to compare the CP violating invariant / = Im tr Mi • • • M n with the "size" (in some 
sense) of the mass matrices involved in its definition. The bare quantity I cannot give this 
measure because with a rescaling M, — > AMj it transforms as I — > X n I (because it has 
dimensions [I] = m n ) and in absence of anomalies, physics should not change with such a 
scale transformation. Thus we have to divide / by a quantity Af with dimensions [AT] = m n , 
invariant under a change of weak quark basis and which in some sense measures the "size" of 
the mass matrices M\ , . . . , M n . The only nonsingular invariant measure of the size of a matrix 
A is its norm \ \A\ \ = (tr j^t) 1 / 2 (the other invariants of A, for instance its determinant, can 
be zero without A being identically zero and can lead to singularities in the definition). The 
simplest choice is then Eq. (^). 

It is apparent that there is no set of parameters (masses, mixings, phases) that maximizes 
all the quantities /. So, any definition of maximal CP violation based on reduced invariants 
will have the arbitrariness in the choice of the reduced invariants that are maximized. The 
most obvious possibility is to choose the lowest order (lowest number of mass matrices in 
the numerator) reduced invariant, because in general higher order reduced invariants involve 
higher powers of the mixing angles and therefore have smaller maximum values than the 
lowest order ones. Having this in mind, we will find the maximum value of the lowest order 
(nontrivial) reduced invariants for some models. 



I = 



Im tr Mi ■ ■ ■ M, 



(3) 



(tr MiMj) 1 / 2 • • • (tr M n Mly/ 2 ' 



3 



As we have seen before, I\ in Eq. (Q) is the lowest order nontrivial invariant in the SM. 
Written as a function of physical parameters, 



h = -{m 2 c - m 2 u )(mf - m 2 u )(m 2 - m 2 c )(m 2 s - m 2 d ){m 2 b - m 2 d )(m 2 b - m 2 s ) 

x Im V ud V; d V cs V: s . (4) 

and the reduced invariant I\ decouples into a factor involving up quark masses only, another 
one with down quark masses and a third one involving CKM matrix elements: 

j _ (mg - m 2 u ){m 2 - m 2 u )(m 2 - m 2 ) ^ (m 2 - m\){m 2 - m 2 )(m 2 b - m 2 ) 
{rn\ + m% + m 2 ) 3 (m d + m 2 s + m 2 ) 3 

x Im V ud V c * d V cs V: s . (5) 

The maximum value of the third factor in Eq. (|5|) was found in Ref. || to be l/6\/3- It is 
straightforward to show that the maximum value of the up (down) factor occurs when one 
mass is zero and the ratio of the other two is (1 + \^3)/V2, and is (surprisingly) also l/6\/3- 
Then, the maximum value of I\ is (l/6\/3) 3 - 

This result can be compared with the hypothetical case that the top quark did not exist, 
i. e. the bottom quark was a vector-like SU(2)i singlet. In this model, the up quark mass 
matrix M u has dimensions 2x2 and the 3x3 down quark mass matrix A4 d divides in two 
submatrices 

M d is a 2 x 3 matrix connecting the two left-handed SU(2)^ doublets with the three right- 
handed SU(2)i singlets; m d is a 1 x 3 matrix connecting the left-handed SU(2)x singlet with 
the three right-handed SU(2) l singlets. (For a detailed description of the models with quark 
singlets see Ref. [13]). The lowest order nontrivial invariant of this model is [12] 

I[ = Im tr H u H d h d h\ 

= -(m 2 c -m 2 J(m 2 s -m 2 d )(ml-m 2 d )(ml-m 2 s )ImV ud V* d V cs V* s , (7) 



with H x = M X M^, h d = M d m} d . It has the same dependence in down quark masses and 
CKM matrix elements than I\ in Eq. (Q). The corresponding reduced invariant is = 
/{/(tr H u )(ti H d ) 2 (tv m d m d ), and its maximum value is 1/16, roughly 70 times larger than 
the maximum value of I\ in the SM. It is worthwhile to note that a naive analysis involving 
only the rephasing invariant CKM factor in Eqs. (|]J7|) yields in both models a maximum value 
of 1/6 \/3- This does not reflect the fact that in the latter model the effects of CP violation 
can be larger than in the SM. The formulation in terms of invariants is more adequate because 
although in both models the CKM matrix is a 3 x 3 unitary matrix (in the latter case only 
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the upper 2x3 submatrix is observable) with one CP violating phase, the invariants are 
different in each case and the observables depend on the masses and CKM matrix elements 
through the invariants. 

These models with small CP violation effects can be compared with models with left- 
right symmetry, where CP violation can be very important. In these models, even with one 
generation of quarks CP violation is possible. The lowest order nontrivial invariant is fl^| 

/(' = Im tr M U M\ . (8) 

In the one generation case, M u and Md are complex numbers m u and rrid respectively, and 
the reduced invariant is I" = Im (m u rn* d ) /\m u \ \md\, which has a maximum value of 1 when 
the relative phase between m u and is ir/2. In the three generation case /[' also reaches 
the value 1 when the up and down mass matrices are equal up to an overall factor. 

We can go a step further in our analysis and try to define quantities that characterize the 
average size of CP violation in a specific model, independently of the actual size (i. e. when 
one inserts the known values of masses, mixings, etc.) of CP violation. The significance of 
this will be clear later. For this purpose we define a vector space X whose elements g are 
the reduced invariants of a given model, which form a basis in this space, plus real linear 
combinations of them. We can regard the elements of I as functions of k real variables, with 
k twice the total number of (complex) entries in the quark mass matrices of the model. (For 
instance, in the SM k = 2 x (3x3 + 3x3) = 36.) We consider the vector subspace lo C I 
consisting of the elements of X which vanish identically for all values of its k variables (the 
reduced invariants / oc Im tr M with M hermitian are some of the elements of Iq). Then we 
construct the vector space Z' = T/Tq. In this space, but not in X, we can define a norm 

| | "| | =(?k/ s |9(I)l S 4 ' (9) 

where x denotes the k variables of g, S is a subset of R fe which we will specify later and V(S) 
is the volume of S. For the moment, we can assume S to be a ball centered at the origin and 
of arbitrary radius r. With this definition, we see: (i) ||<7|| > 0; (ii) If ||<?|| = 0, g(x) = in S, 
but our definition of S implies g(x) = Vx S R fc and g G Tq; (hi) \\ag\\ = \a\ \ \g\\ with a a 
real number; and (iv) +52II < 115111 + 115211 as a direct consequence of the same property 
for functions of k variables. Thus || • || defines a norm on Z' . Note that the property of the 
reduced invariants — 1 < I < 1 implies that the integrand is bounded, so the integral is well 
defined. At this point it is important to note that if we choose another normalization for the 
invariants, we can define the vector spaces X, Iq, X' in the same way and the operation || • || 
is still a norm although the modulus of the vectors differs. 

Before we give a significance to the norm we will specify the set S. This is best explained 
with an example. Let us consider the SM with quark mass matrices M u = (rriij), = (n^). 
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The elements of the set S discussed above are points x £ R fc with components Re m«, 
Im rriij, Re riij, Im mj, i,j = 1, 2, 3 such that Re 2 m,jj + Im 2 mjj + Re 2 ny + Im 2 ?^ < r 2 . 
This set is not invariant under weak basis transformations, so it is more convenient to define 
S = Re 2 mjj + lm 2 mij < r 2 , Re 2 njj + Im 2 njj < r 2 } i.e. the set S is the direct 

product of two 18-dimensional balls centered at the origin. The norm \\g\\ of a vector is 
independent of the value of r because g is dimensionless. In a general model, the set S can 
be taken to be the direct product of a ball of radius 1 centered at the origin containing the 
parameters for each mass matrix. Then, the significance of ||/|| is: (i) Choose a point in the 
parameter space of the theory; (ii) Calculate the measure of CP violation in this theory with 
these parameters, as given by /; (iii) Average over a symmetric set of points containing all 
possible values of the parameters. The result is an "average measure" of CP violation in the 
model, as given by /. The interest of these average measures constructed as norms of reduced 
invariants is that they tell us a priori and without knowledge of the values of the parameters 
in the theory whether CP violation effects are expected to be large or small. Of course, the 
size of CP violation effects is given by the actual parameters, but these average measures 
give us the global behavior of the theory. Comparing the norms of the reduced invariants of 
two theories we can state in which theory CP violation is more important. The numerical 
calculation of the norms of I±, I[, I'{ (the latter for three generations of quarks) yields 

||Ji|| = 1.63 -1(T 5 , 
= 2.51 -10 -4 , 

= 0.4003. (10) 

The evaluation of the multidimensional integral of Eq. (|9|) to obtain Eqs. (|io| ) is a difficult 
task because the integrand is very oscillatory (specially in the first case) and we integrate 
over a huge number of dimensions (36, 26 and 36 respectively), so we have used a Monte 
Carlo method and the errors are on the last decimal place. We see that ||/i|| < \\I'i\\ and 
ll^ill < W\\\ as expected from the hierarchy of their maximum values. However, while in 
the model with left-right symmetry we have | \I'{ \ ) m ax — 0.4, in the other two models 
we have ||ii||/(Ii) max ~ 1.8 • 10~ 2 , 1 \/(I[) max ~ 4.0 • 10~ 3 (so the two first integrals are 
strongly peaked and difficult to calculate). This means that the average size of CP violation 
in the SM and the model with vector-like bottom is much less than the maximum allowed 
by the model, i. e. by the quark content and symmetries of the Lagrangian. These results 
are expected to be general and not sensitive to the choice of normalization. Within the SM, 
using the experimental values of masses and mixings we find h < 9.1 • 10~ 13 . Thus, in the 
SM not only CP violation is much less than maximal, but is also much less than the average. 
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